Locating Boosted Kerr and Schwarzschild Apparent Horizons

We describe a finite-difference method for locating apparent horizons and illustrate its capabilities on boosted Kerr and Schwarzschild black holes. Our model spacetime is given by the Kerr-Schild metric. We apply a Lorentz boost to this spacetime metric and then carry out a 3+1 decomposition. The result is a slicing of Kerr/Schwarzschild in which the black hole is propagated and Lorentz contracted. We show that our method can locate distorted apparent horizons efficiently and accurately.Comment: Submitted to Physical Review D. 12 pages and 22 figure

Extended Lifetime in Computational Evolution of Isolated Black Holes

Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint equations remain solved under the action of the evolution, and one approach is to simply monitor them ({\it unconstrained} evolution). The problem of the 3-d computational simulation of even a single isolated vacuum black hole has proven to be remarkably difficult. Recently, we have become aware of two publications that describe very long term evolution, at least for single isolated black holes. An essential feature in each of these results is {\it constraint subtraction}. Additionally, each of these approaches is based on what we call "modern," hyperbolic formulations of the Einstein equations. It is generally assumed, based on computational experience, that the use of such modern formulations is essential for long-term black hole stability. We report here on comparable lifetime results based on the much simpler ("traditional") $\dot g$ - $\dot K$ formulation. We have also carried out a series of {\it constrained} 3-d evolutions of single isolated black holes. We find that constraint solution can produce substantially stabilized long-term single hole evolutions. However, we have found that for large domains, neither constraint-subtracted nor constrained $\dot g$ - $\dot K$ evolutions carried out in Cartesian coordinates admit arbitrarily long-lived simulations. The failure appears to arise from features at the inner excision boundary; the behavior does generally improve with resolution.Comment: 20 pages, 6 figure

Initial Data and Coordinates for Multiple Black Hole Systems

We present here an alternative approach to data setting for spacetimes with multiple moving black holes generalizing the Kerr-Schild form for rotating or non-rotating single black holes to multiple moving holes. Because this scheme preserves the Kerr-Schild form near the holes, it selects out the behaviour of null rays near the holes, may simplify horizon tracking, and may prove useful in computational applications. For computational evolution, a discussion of coordinates (lapse function and shift vector) is given which preserves some of the properties of the single-hole Kerr-Schild form

Generic Tracking of Multiple Apparent Horizons with Level Flow

We report the development of the first apparent horizon locator capable of finding multiple apparent horizons in a ``generic'' numerical black hole spacetime. We use a level-flow method which, starting from a single arbitrary initial trial surface, can undergo topology changes as it flows towards disjoint apparent horizons if they are present. The level flow method has two advantages: 1) The solution is independent of changes in the initial guess and 2) The solution can have multiple components. We illustrate our method of locating apparent horizons by tracking horizon components in a short Kerr-Schild binary black hole grazing collision.Comment: 13 pages including figures, submitted to Phys Rev

Three-dimensional adaptive evolution of gravitational waves in numerical relativity

Adaptive techniques are crucial for successful numerical modeling of gravitational waves from astrophysical sources such as coalescing compact binaries, since the radiation typically has wavelengths much larger than the scale of the sources. We have carried out an important step toward this goal, the evolution of weak gravitational waves using adaptive mesh refinement in the Einstein equations. The 2-level adaptive simulation is compared with unigrid runs at coarse and fine resolution, and is shown to track closely the features of the fine grid run.Comment: REVTeX, 7 pages, including three figures; submitted to Physical Review